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METHOD AND SYSTEM FOR PREDICTIVE PHYSIOLOGICAL GATING OF 

RADIATION THERAPY 

BACKGROUND OF THE INVENTION 
Field Of The Invention 

The present invention relates to medical methods and systems. More particularly, 
the invention relates to a method and system for physiological gating of radiation therapy. 

Background 

Radiation therapy involves medical procedures that selectively expose certain areas 
of a human body, such as cancerous tumors, to high doses of radiation. The intent of the 
radiation therapy is to irradiate the targeted biological tissue such that the harmful tissue is 
destroyed. In certain types of radiotherapy, the irradiation volume can be restricted to the 
size and shape of the tumor or targeted tissue region to avoid inflicting unnecessary 
radiation damage to healthy tissue. For example, conformal therapy is a radiotherapy 
technique that is often employed to optimize dose distribution by conforming the treatment 
volume more closely to the targeted tumor. 

Normal physiological movement represents a limitation in the clmical planning and 
delivery of conventional radiotherapy and conformal therapy. Normal physiological 
movement, such as respiration or heart movement, can cause a positional movement of the 
tumor or tissue region undergoing irradiation. If the radiation beam has been shaped to 
conform the treatment volume to the exact dimensions of a tumor, then movement of that 
tumor during treatment could result in the radiation beam not being sufficiently sized or 
shaped to fully cover the targeted tumoral tissue. 

One approach to this problem mvolves physiological gating of the radiation beam 
during treatment, with the gating signal synchronized to the movement of the patient's 
body. In this approach, instruments are utilized to measure the physiological state and/or 
movement of the patient. For example, respiration has been shown to cause movements in 
the position of a lung tumor in a patient's body. If radiotherapy is being applied to the limg 
tumor, then a temperature sensor, strain gauge, preumotactrograph, or optical imaging 
system can be utilized to measure the patient's respiration cycle. These instruments can 
produce a signal indicative of the movement of the patient during the respiratory cycle. The 
radiation beam can be gated based upon certain threshold amplitude levels of the measured 
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respiratory signal, such that the radiation beam is disengaged or stopped during particular 
time points in the respiration signal that correspond to excessive movement of the lung 
tumor. 

Known approaches to physiological gating of radiation therapy are reactive, that is, 
5 known approaches utilize gating methods that slavishly react to measured levels of 

physiological movements. One drawback to reactive gating systems is that the measured 
physiological movement may involve motion that that is relatively fast when compared to 
the effectively operating speeds of gating system components. Thus, a purely reactive 
gating system may not be able to react fast enough to effectively gate the applied radiation. 
10 For example, the gating system may include a switch for gating the radiation treatment, in 
which the switch requires a given time period At to fully engage. If the switching time 
period At is relatively slow compared to the measured physiological motion cycle, then a 
system employmg such a switch in a reactive manner may not be able to effectively gate the 
application of radiation at appropriate tune points during the radiation therapy. 
1 5 Therefore, there is a need for a system and method to address these and other 

problems of the related art. There is a need for a method and system for physiological 
gating which is not purely reactive to measure physiological movement signals. 
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SUMMARY OF THE INVENTION 
The present invention provides an improved method and system for physiological 
gatmg for radiation therapy. According to an aspect, the invention comprises a method and 
system for detecting and predictably estimating regular cycles of physiological activity or 
5 movements. Another aspect of the invention is directed to predictive actuation of gating 
system components. Yet another aspect of the invention is directed to physiological gating 
of radiation treatment based upon the phase of the physiological activity. 

These and other aspects, objects, and advantages of the invention are described 
below in the detailed description, drawings, and claims. 

10 

BRIEF DESCRIPTION OF THE DRAWINGS 
The accompanying drawings are included to provide a further understanding of the 
invention and, together with the Detailed Description, serve to explain the principles of the 
invention. 

15 Fig. 1 depicts the components of a system for physiological gatmg according to an 

embodiment of the invention. 

Fig. 2 depicts an example of a respiratory motion signal chart. 

Fig. 3 depicts a motion signal chart and a gating signal chart. 

Fig. 4 is a flowchart showing process actions performed in an embodiment of the 
20 invention. 

Fig. 5 is a flowchart showing process actions for detecting and predicting estimation 
of regular physiological movements. 

Fig. 6a depicts a side view an embodiment of a camera that can be utilized in the 
invention. 

25 Fig. 6b depicts a front view of the camera of Fig. 6a. 

Fig. 7a depicts a retro-reflective marker according to an embodiment of the 
invention. 

Fig. 7b depicts a cross-sectional view of the retro-reflective marker of Fig. 7a. 
Fig. 8 depicts an apparatus for making a retro-reflective marker. 
30 Fig. 9 depicts a phase chart synchronized with a gating signal chart. 

Fig. 10 depicts an embodiment of a cylindrical marker block. 
Fig. 1 1 depicts an embodiment of a hemispherical marker block. 
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Fig. 12 is a diagram of a computer hardware system with which the present 
invention can be implemented. 

DETAILED DESCRIPTION 

An aspect of the present invention comprises a method for detecting and predictively 
estimating regular cycles of physiological activity or movement. The method of the 
invention can be employed for any regular physiological activity, including for example, the 
respiratory or cardiac cycles. 

In operation, one or more sets of data representative of the physiological activity of 
interest are collected for the patient. For example, an electrocardiograph can be employed 
to generate data representative of the cardiac cycle. To generate data representative of the 
respiration cycle, a temperature sensor, strain gauge or preumotactrograph can be employed. 
Other instruments or mechanisms can be employed within the scope of the invention to 
obtain sets of data representative of physiological activity or movements. 

Fig. 1 depicts the components of an embodiment of a system 100 for physiological 
gating of radiation therapy, in vAdch data representative of physiological activity is 
collected with an optical imaging apparatus. System 100 comprises a radiation beam source 
102 (such as a conventional linear accelerator) which is positionally configured to direct a 
radiation beam at a patient 106 located on treatment table 104. A switch 1 16 is operatively 
coupled to the radiation beam source 102, Switch 1 16 can be operated to suspend the 
application of the radiation beam at patient 106. In an embodiment, switch 1 16 is part of 
the mechanical and electrical structure of radiation beam source 102. Alternatively, switch 
1 16 comprises an external apparatus that is connected to the control electronics of radiation 
beam source 102. 

An optical or video image apparatus, such as video camera 108, is aimed such that 
as least part of the patient 106 is within the camera's field of view. Camera 108 monitors 
patient 106 for motion relating to the particular physiological activity being measured. For 
example, if respiration movements of the patient is being monitored, then camera 108 can 
be configured to monitor the motion of the patient's chest. According to an embodiment, 
camera 108 is placed with its axis at approximately 45 degrees to the longitudinal axis of 
the patient 106. For measurement of respiration activity that could result in about 3-5 mm 
of chest motion, the video image field of view is preferably set to view an approximately 20 
cm by 20 cm area of the patient's chest For purposes of illustration only, a smgle camera 
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108 is shown in Fig. 1. However, the number of cameras 108 employed in the present 
invention can exceed that number, and the exact number to be used in the invention depends 
upon the particular application to which it is directed. 

In an embodiment, one or more illumination sources (which are infrared sources in 
5 the preferred embodiment) project light at the patient 1 06 oh treatment table 1 04. The 
generated light is reflected from one or more landmarks on the patient's body. The camera 
108, which is directed at patient 106, captures and detects the reflected light from the one or 
more landmarks. The landmarks are selected based upon the physiological activity being 
studied. For example, for respiration measurements, landmarks are selected from one or 
1 0 more locations on the patient's chest. 

The output signals of camera 108 are sent to a computer 1 10 or other type of 
processing unit having the capability to receive video images. According to a particular 
embodiment, computer 1 10 comprises an Intel Pentium-based processor running Microsoft 
Windows NT and includes a video frame grabber card having a separate channel for each 
1 5 video source utilized in the system. The images recorded by camera 1 08 are sent to 

computer 1 10 for processing. If camera 108 produces an analog output, the frame grabber 
converts the camera signals to a digital signal prior to processing by computer 1 10. Based 
upon the video signals received by computer 1 10, control signals can be sent from computer 
1 1 0 to operate switch 116. 
20 According to one embodiment, one or more passive markers 1 14 are located on the 

patient in the area to be detected for movement. Each marker 1 14 preferably comprises a 
reflective or retro-reflective material that can reflect light, whether in the visible or invisible 
wavelengths. If the illumination source is co-located with camera 108, then marker 1 14 
preferably comprises a retro-reflective material that reflects light mostly in the direction of 
25 the illumination source. Alternatively, each marker 1 14 comprises its own light source. 
The marker 1 14 is used in place of or in conjunction with physical landmarks on the 
patient's body that is imaged by the camera 108 to detect patient movement. Markers 1 14 
arc preferably used instead of body landmarks because such markers 1 14 are easier to detect 
and track via the video image generated by camera 108. Because of the reflective or retro- 
30 reflective qualities of the preferred markers 1 14, the markers 1 14 inherently provide greater 
contrast in a video unage to a light detecting apparatus such as camera 108, particularly 
when the camera 108 and illumination source are co-located. 



5 



Docket No. 98-43 



Utilizing a video or optical based system to track patient movement provides several 
advantages. First, a video or optical based system provides a reliable mechanism for 
repeating measurement results between uses on a given patient. Second, the method of the 
invention is noninvasive, and even if markers are used, no cables or connections must be 
5 made to the patient. Moreover, if the use of markers is impractical, the system can still be 
utilized without markers by performing measurements of physiological activity keyed to 
selected body landmarks. Finally, the method of the invention is more accurate because it is 
based upon absolute measurement of external anatomical physical movement. 

A possible inefficiency in tracking the markers 1 14 is that the marker may appear 
1 0 anywhere on the video frame, and all of the image elements of the video frame may have to 
be examined to determine the location of the marker 1 14. Thus, in an embodiment, the 
initial determination of locations for the marker 1 14 involves an examination of all of the 
image elements in the video frame. If the video frame comprise 640 by 480 image 
elements, then all 307200 (640'*'480) image elements are initially examined to find the 
1 5 location of the markers 1 14. 

For real-time tracking of the marker 114, examining every image element for every 
video frame to determine the location of the marker 1 14 in real-time could consume a 
significant amount of system resources. Thus, in an embodiment, the real-time tracking of 
marker 1 14 can be facilitated by processing a small region of the video frame, referred to 
20 herein as a "tracking gate", that is placed based on estimation of the location of the already- 
identified marker 1 14 in the video frame. The previously determined location of a meirker 
114 defined in the previous video frame is used to define an initial search range (i.e., the 
tracking gate) for that same marker in real-time. The tracking gate is a relatively small 
portion of the video frame that is centered at the previous location of the marker 1 14. The 
25 tracking gate is expanded only if it does not contain the new location of the marker 1 14. As 
an example, consider the situation when the previously determined location of a particular 
marker is image element (50,50) in a video frame. If the tracking gate is lunited to a SO by 
50 area of the video frame, then the tracking gate for this example would comprise the 
image elements bound within the area defined by the coordinates (25,50), (75,50), (50,25), 
30 and (50,75). The other portions of the video frame are searched only if the marker 106 is 
not found within this tracking gate. 

The video image signals sent from camera 108 to computer 1 10 are used to generate 
and track motion signals representative of the movement of marker 1 14 and/or landmark 
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structures on the patient's body. Fig. 2 depicts an example of a motion signal chart 200 for 
respiratory movement that contains information regarding the movement of marker 1 14 
during a given measurement period. The horizontal axis represents points in time and the 
vertical axis represents the relative location or movement of the marker 1 14. According to 
an embodiment, the illustrated signal in Fig. 2 comprises a plurality of discrete data points 
plotted along the motion signal chart 200. 

An important aspect of physiological gating of radiotherapy is the determination of 
the boundaries of the "treatment intervals" for applying radiation. For gating purposes, 
threshold points can be defined over the amplitude range of the motion signal to determine 
the boundaries of the treatment mtervals. Motion of the patient that fall outside the 
boundaries of the treatment intervals correspond.to movement that is predicted to cause 
unacceptable levels of movement to the tumor or tissue targeted for irradiation. According 
to an embodiment, the treatment intervals correspond to the portion of the physiological 
cycle in which motion of the clinical target volume is minimized. Other factors for 
determining the boundaries of the treatment intervals include identifying the portion of the 
motion signals involvmg the least movement of the target volume or the portion of the 
motion signal involving the largest separation of the target volume from organs at risk. 
Thus, the radiation beam pattern can be shaped with the minimum possible margin to 
account for patient movement. 

Radiation is applied to the patient only when the motion signal is within the 
designated treatment intervals. Referring to Fig. 3, depicted are examples of treatment 
mtervals, indicated by signal range 302, that has been defined over the motion data shown 
in motion signal chart 200. In the example of Fig. 3, any movement of the measured body 
location that exceeds the value of 0.8 (shown by upper boundary line 304) or which moves 
below the value of 0.0 (shown by lower boundary line 306) falls outside the boundaries of 
the treatment intervals. 

Shown in Fig. 3 is an example of a gating signaLchart 300 that is aligned with 
motion signal chart 200. Any motion signal that falls outside the treatment interval signal 
range 302 results in a "beam hold" gating signal threshold 310 that stops the application of 
radiation to the patient. Any motion signal that is within the treatment interval signal range 
302 results in a "beam on" gating signal threshold 312 that allows radiation to be applied to 
the patient. In an embodiment, digital signals that represent the information shown in 
motion signal chart 200 are processed by computer 1 10 and compared to the threshold 
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levels of the treatment interval signal range 302 to generate gating signal thresholds 3 10 and 
312. Alternatively, gating signal thresholds 310 and 312 can be obtained by feeding analog 
motion signals to a comparator to be compared with analog threshold signals that 
correspond to treatment interval signal range 302. In any case, gating signal thresholds 310 
and 3 12 are generated by computer 110 and are applied to the switch 1 16 that controls the 
operation of radiation beam source 102 (Fig. 1) to stop or start the application of a radiation 
beam at patient 1 06. 

Fig. 4 is a flowchart of the process actions performed in an embodiment of the 
invention. The first process action is to define boundaries for the treatment intervals over 
the range of motion signals to be detected by a camera (402). As indicated above, any 
motion that fall outside the boundaries of the treatment intervals correspond to motion that 
is predicted to result in xmacceptable levels of movement of the tumor or tissue targeted for 
irradiation. An optical or video imaging system, such as a video camera, is used to 
measure the physiological motion of the patient (404), and the output signals of the optical 
or video imaging system are processed to compare the measured motion signals with the 
threshold boundaries of the treatment intervals (406). 

If the motion signal is outside the boundaries of the treatment intervals, then a 
"beam off gating signal threshold is applied to a switch that is operatively coupled to the 
radiation beam source (408). If the radiation beam source is presently irradiating the patient 
(410), then the switch setting is operated to hold or stop the radiation beam (41 1), The 
process then returns back to process action 406. 

If the motion signal is within the boxmdaries of the treatment intervals, then a "beam 
on" gating signal threshold is produced (412) and is applied to a switch that is operatively 
coupled to the radiation beam source. If the radiation beam source is presentiy not being 
applied to the patient (413), then the switch setting is operated to turn on or apply the 
radiation beam source to irradiate the patient (414). The process then returns back to 
process action 406. 

According to one embodiment, the radiation beam source can be disengaged if a 
significant deviation is detected in the regular physiological movements of tiie patient. 
Such deviations can result fi-om sudden movement or coughing by the patient The position 
and/or orientation of the targeted tissue may unacceptably shift as a result of this deviation, 
even though the amplitude range of the motion signal still falls within the boundaries of the 
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treatment intervals during this deviation. Thus, detection of such deviations helps define the 
appropriate time periods to gate the radiation treatment. 

The present invention provides a method for detecting and predictively esthnating a 
period of a physiological activity. In effect, the present invention can "phase lock" to the 
5 physiological movement of the patient. Since the gating system phase locks to the 

movement period, deviations from that period can be identified and appropriately addressed. 
For example, when gating to the respiratory cycle, sudden movement or coughing by the 
patient can result in deviation fi-om the detected period of the respiration cycle. The 
radiation treatment can be gated during these deviations from the regular period. The 
10 present invention also provides a method for predictively estimating the period of the 
subsequent physiological movement to follow. . 

Fig. 5 is a process flowchart of an embodiment of the invention to perform 
predictive estimation and detection of regular physiological movement cycles. In process 
action 502, an instrument or system (such as system 100 from Fig. 1) is employed to 
15 generate data signals representative of the physiological activity of interest. In an 

embodiment, the data signals comprises a stream of digital data samples that collectively 
form a signal wave pattern representative of the physiological movement \mder 
examination. A number of discrete measurement samples are taken for the physiological 
activity during a given time period. For example, in an embodiment of the invention 
20 directed towards respiratory measurement, approximately 200-210 data samples are 
measured for each approximately 7 second time interval. 

In process action 504, pattern matching analysis is performed against the measured 
data samples. In an embodiment, the niost recent set of data samples for the physiological 
signal is correlated against an immediately preceding set of data samples to determine the 
25 period and repetitiveness of the signal. An autocorrelation fimction can be employed to 
perform this pattern matching. For each new sample point of the physiological motion or 
physiological monitoring sensor signal, the process computes the autocorrelation fimction of 
the last n samples of the signal, where n corresponds to approximately 1.5 to 2 signal 
breathing periods. The secondary peak of the autocorrelation fimction is then identified to 
30 determine the period and repetitiveness of the signal. 

In an alternate embodiment, an absolute difference fimction is used instead of an 
autocorrelation fimction. Instead of secondary peak, a secondary minimum in the absolute 
difference is searched for. For each hew sample point of the physiological motion or 
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physiological monitoring sensor signal, the process computes the minimxim absolute 
difference between fte two sets of data over a range of overlapping data samples. The 
secondary minimum corresponds to the data position that best matches the recent set of data 
samples with the preceding set of data samples. 

Yet another altemate embodiment performs a pattern matching based upon a model 
of the physiological activity being measured. The model is a dynamic representation of the 
physiological motion or physiological monitoring sensor signal for that physiological 
activity. The latest set of data samples is matched against the model to estimate parameters 
of the repetitive process. 

Pattern matching xising the measured physiological signal (504) provides 
information regarding the degree of match, as well as a location of best match for the 
repetitive process. If an autocorrelation function is employed in process action 504, then 
the relative strength of secondary peak provides a measure of how repetitive the signal is. A 
threshold range value is defined to provide indication of the degree of match between the 
two sets of data samples. If the strength of the secondary peak is within the defmed 
threshold range (process action 508), then the degree of match indicates that the signal is 
repetitive, and the secondary peak location provides an estimate of the signal period. If an 
absolute difference function is used in process action 504, then the relative value of the 
secondary minimum provides a measure of how repetitive the signal is. If the value of the 
secondary minimum meets a defined threshold range (508), then the degree of match 
indicates that the signal is repetitive, and the secondary minimum location provides an 
estimate of the signal period. 

If the correlation value of the secondary peak or secondary minimmn does not meet 
the defined threshold range, then a deviation from the regular physiological activity is 
detected, thereby indicating an irregularity in the regular physiological movement of the 
patient (510). This irregularity could result, for example, fi:om sudden movement or 
coughing of the patient. In an embodiment, this detected irregularity results in the 
generation of a "beam hold" signal that suspends the application of radiation at the patient. 

If the degree of match indicates repetitiveness, the point of best match is tested to 
determine if the period is within a reasonable range. The location of the secondary peak or 
secondary minimum provides an estimate of the period of the physiological activity. In an 
embodiment, the point of best match is compared to a threshold range (509). If the point of 
best match does not fall within the threshold range, than a deviation firom regular 
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physiological activity is detected (510). If the point of best match falls within the threshold 
range, then the signal is accepted as being repetitive (512). 

The estimate of tiie period based on the point of best match can be used to predict 
the period and waveform parameters of the next set of data samples for the signal (514), 
Note that process actions 504, 508, and 509 test for repetitiveness based upon a plurality of 
data samples over a range of such samples. However, in some circumstances, a significant 
deviation fi-om normal physiological movement may actually occur within the new or most 
recent data sample(s) being analyzed, but because the overall set of data samples indicates 
repetitiveness (e.g., because of averaging of absolute differences over the range of data 
samples being compared), process actions 504, 508, and 509 may not detect the deviation. 
To perform a test for rapid deviation, the predicted value from process action 514 is 
compared with the next corresponding data sample (515), If the predicted value does not 
match the actual data sample value witiiin a defined threshold range, then deviation is 
detected (510). If a comparison of the predicted and actual data sample values fall within 
the defined threshold range, then repetitiveness is confirmed, and deviation is not detected 
for that data sample range (5 1 6). 

In an embodiment, the first time the process of Fig. 5 is performed, the pattern 
matching process action (504) is performed over the entire range of data samples. 
Thereafter, the pattern matching process action can be performed over a limited search 
interval, which is defined by the results of the prior immediate execution of the process. 
For example, the predicted value from process action 514 can be used to define the location 
of the search interval for the next set of data samples. However, if process action 508, 509, 
and 514 detect deviation based upon analysis of the initial search interval, then the search 
interval can be expanded to ensure that a deviation has actually occurred. The process of 
Fig. 5 can be repeated with the increased search interval to attempt to find a point of best 
match outside of the initial search interval. In an embodiment, this increased search interval 
comprises the entire range of data samples. Alternatively, the increased search interval 
comprises only an expanded portion of the entire range of data samples. 

According to an embodiment of the invention, physiological gating can be 
performed based upon the phase of the physiological activity, rather than its amplitude. 
This is in contrast to the example shown in Fig. 3, in which the amplitude of the 
physiological movement signal defines the boundaries of the treatment intervals for the 
application of radiation. 
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Referring to Fig. 9, depicted is an example of a chart 900 showing the phase 

progression of a physiological movement signal. Treatment interval range 902 have been 

defined over the phase chart 900. In the example of Fig. 9, the boundaries of the treatment 

interval range 902 are defined by the phase of the detected signal. Radiation is applied to 

the patient only when the phase of the physiological movement signal falls within the 

boundaries of the treatment interval range 902. Fig. 9 depicts examples of treatment 

interval range 902 having boundaries that span from 30 degrees to 300 degrees. Thus, the 

applied radiation to the patient is suspended or stopped when the phase of the physiological 

movement signal is between 301 degrees to 29 degrees. 

Shown in Fig. 9 is a gating signal chart 906 that is aligned with phase chart 900. A 

"beam hold" signal threshold 910 results if the phase of the physiological movement signal 

falls outside the treatment interval range 902. A "beam on'' signal threshold 912 results if 

the phase of the physiological movement signal falls within the boundaries of the treatment 

interval range 902. The "beam on" and "beam hold" signal thresholds 910 and 912 are 

applied to a switch 1 16 that operatively controls the operation of a radiation beam source 

102 (Fig. 1). If radiation is being applied to the patient, application of the "beam hold" 

signal threshold 910 triggers switch 1 16 to suspend or stop the application of radiation. If 

radiation to the patient is not being applied, application of the "beam on" signal threshold 

912 triggers the application of radiation to the patient. 

The predictive qualities of the present invention permits operation of a gating system 

even if the measured physiological movement involves motion that that is relatively fast 

when compared to the effectively operating speeds of gating system components. As just 

one example, consider a gating system that includes a switch for gating the radiation 

treatment, in which the switch requires a known time period At to fully engage. If the 

switching time period At is relatively slow compared to the measured physiological motion 

cycle, then a system employing such a switch in a reactive manner may not be able to 

effectively gate the application of radiation at the patient. 

The present invention allows predictive triggering of switch 1 16 to compensate for 

the amount of time At required to fiilly engage the switch. A predicted period for a 

physiological activity can be obtained by employing the process of Fig. 5. A treatment 

interval range is defined over a portion of the period of the physiological activity. Based 

upon the time At required to fully actuate the switch 1 16, the switch 1 16 can be pre-actuated 

by this time period At prior to the time of the boundary of the treatment interval, so that the 

12 



Docket No. 



time for full actuation of the switoh 116 coincides with the boundary of the treaUnent 
interval. Thus, the radiation can be effectively gated at the boundaries of the treatment 
interval, regardless of the operating speeds of the switch 116. The same procedure can also 
be employed to compensate for the operating speeds of other gating system components. 
5 The following is an embodiment of the invention coded in the Visual Basic 

programming language. The following program code is directed to a process for detecting 
and predictively estimating the respiration cycle period using the absolute difference 
function: 

1 0 Public Function Predict(By Val i As Long, By Val Range As Long, Period As Double, 
MinAbsDiffAs Double, Diff As Double) As Double 

Dim j As Long, StartJ As Long, CurrJ As Long 
Dim k As Long, MaxK As Long 
1 S Dim AbsDiff As Double 

Dim NorniAbsDiff As Double, n As Long 
k = Period - Range 
MinAbsDifif = 10000000# 

StartJ = TimeRefIdxBuf((i - 201 + BufLength) Mod BufiLength) 
20 CurrJ = TuneRefIdxBuf((i - 1 + BufLength) Mod BufLength) 
Do 
j = StartJ 
AbsDifF= 0# 
n = 0 
25 Do 

AbsDiff = AbsDifF+ Abs(SigBuf(SigRefldxBuf(i)) - SigBufl[SigRefIdxBufi[a + k + 
ChartWidth) Mod ChartWidth))) 
n = n+ 1 

j = 0 + 10) Mod ChartWidth 
30 Loop While n <= (200 - k) / 10 

NormAbsDiff = 100 * AbsDiff / (n * Abs(MaxSignal - MinSignal)) 
If NormAbsDiflf <= MinAbsDifif Then 
MinAbsDiff = NormAbsDifif 
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MaxK = k 

End If 

k = k+l 
Loop While k <= Period + Range 
If MaxK >= 40 And MaxK <- 1 50 Then Period = MaxK 

Predict = SigBuf(SigRefldxBuf((CurrJ - Period + ChartWidth) Mod ChartWidth)) 

Diff = 100 * Abs(SigBuf(SigRefldxBuf(CurrJ)) - Predict) / Abs(MaxSignal - MinSignal) 

If MinAbsDiff <= 20 Then 

ProgressBarl. Value = MinAbsDiff 
Else 

ProgressBarl .Value = 20 
End If 

End Function 

In this program code, the variable "i" represents a coxmter or index to the data 
sample being processed. The variable "Range" represents the search range that is to be 
analyzed. If the period of the physiological cycle has already been determined (i.e., from a 
prior execution of this program code), then the variable "Period" comprises the detected 
period. If the period has not yet been determined, then the variable "Period" is set to a 
default value representative of a normal respiration cycle (e.g., the number of data points in 
a normal breathing cycle, which is approximately 95 data samples in an embodiment of the 
invention where approximately 200-210 data samples are obtained over an approximately 7 
second time period). The variable "MinAbsDiff' is the minimum absolute difference value 
over the search range. The variable "Diff * represents a comparison between the actual 
value of a next data sample and the expected value of that next data sample. 

The variables "j", "StartP*, and "CurrJ" are counters or indexes into the data samples 
being processed. The variable "k" is a counter to the search range. The variable "MaxK" 
represents the position in the search range having the minimum absolute difference value. 
The variable "AbsDifif" maintains the sum of the absolute difference values for overlapping 
data samples. The variable **NormaAbsDiff" is the average absolute difference value for a 
particular position in the search range, which is represented as a percentage value. The 
variable "n" is used to track the position of the data samples relative to the search range 
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which is represented as a percentage value. "Predict" is the predicted value that is returned 
by this program code. 

The variable "MinAbsDiff ' is initialized to a high value so that so that any 
subsequent absolute difference value will be smaller than the initialized value. In an 
embodiment, the set of data samples being processed comprises 200 data points. Thus, in 
this program code, the variable "StartJ" is initialized back 201 data samples. The variable 
"CurrJ" is initialized back one data sample. Because a circular array is being used, the 
"BufLength" variable is referenced during the initialization of both "StartJ" and "CurrJ". 

The outer Do loop moves the current and preceding sets of data samples relative to 
each other. The outer Do loop is active while the variable "k" indicates that the program 
code is processing within the search range. In an embodiment, the search range is initially 
set at three positions to either side of a predicted position. The predicted position is based 
upon the period obtained for an immediately prior execution of the program code. If the 
program code has not been executed immediately prior, then a default period value is used. 
If an acceptable minimum absolute difference value is not found within this search range, 
then the search range can be expanded to, for example, 50 positions to either side of the 
predicted position. 

The variable "j" is mitialized to the "StartP' value. The variables "AbsDifiF' and "n" 
are also initialized prior to execution of the inner Do loop. 

The inner Do loop performs the computation of the absolute difference values 
between the present set and prior set of data samples. The variable "AbsDifF' maintains the 
sum of the absolute difference of values for overlapping data samples being compared. 
Note that the number of data samples being analyzed to determine the absolute difference 
values varies based upon the position in the search range being processed. This results 
because different positions in the search range has different numbers of data samples that 
overlap with the previous set of data samples being compared. In the embodiment of this 
program code, the absolute difference function is computed using every 10^ signal sample 
pomt, i.e., a subsampled subtraction is used. Because a circular array is being used, the 
"Chartwidth" variable is referenced during the calculation of "AbsDifP\ 

The variables *'MaxSignal" and "MinSignal" indicate a maximum and minimum 
range for signal values that have previously been sampled. These values can be established, 
for example, during a learning period for the system ui which data samples are obtained for 
a plurality of respiratory cycles. The •*NormAbsDiff' variable holds the average absolute 
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difference value represented as a percentage value based upon the "MaxSignal" and 
"MinSignal" values. 

If the *TNlormAbsDiff'* value is less than or equal to a previously established 
"MinAbsDiflF' value, then the "MinAbsDifT' variable is set to the ''NomiAbsDifP' value. 
The "MaxK" variable is set to the value of "k" if the "MinAbsDifif" value is reset. The 
variable "k" is then incremented, and if the "k" value is still within the search range, then 
the program code returns back to the beginning of the outer Do loop. 

The result of this program code is a candidate minimum absolute difference value 
and a candidate position for the minimum absolute difference. The MaxK value is 
compared to pre-defined threshold values to ensure that it falls within a correct range of 
values for the physiological activity being processed. Thus, in an embodiment, the MaxK 
value is tested to make sure that it is greater than or equal to 40 and less than or equal to 
150 If the mark value meets the threshold range, then the variable "Period" is set to the 
"MaxK" value. The variable "Predict" returns the predicted value for the next set of data 
samples to be processed. The variable "DifF* indicates the difference value between the 
cxirrent data sample value and the predicted data sample value, and is represented as a 
percentage to the "MaxSignal" and "MixSignal" values. 

In an embodiment, an image of a progress bar can be displayed to visually indicate 
the periodicity of the signal samples. Accordmg to the program code, if the "MinAbsDifif" 
value is less than or equal to a 20% difference, then the visual progress bar is updated with 
the computed "MinAbsDiff' value. Otherwise, the visual progress bar displays all other 
"MinAbsDiff" values that exceed a 20% difference as a default value of "20". 

Figs. 6a and 6b depict an embodiment of a camem 108 that can used in the present 
invention to optically or visually collect data representative of physiological movement. 
Camera 108 is a charge-couple device ("CCD") camera having one or more photoelectric 
cathodes and one or more CCD devices. A CCD device is a semiconductor device that can 
store charge in local areas, and upon appropriate control signals, transfers that charge to a 
readout point When light photons from the scene to be images are focussed on the 
photoelectric cathodes, electrons are liberated in proportion to light intensity received at the 
camera. The elections are captured in charge buckets located within the CCD device. The 
distribution of captured electrons in the charge buckets represents the unage received at the 
camera. The CCD transfers these electrons to an analog-to-digital converter. The output of 
the analog-to-digital converter is sent to computer 410 to process the video image and to 
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calculate the positions of the retro-reflective markers 406. According to an embodiment of 
the invention, camera 108 is a monochrome CCD camera having RS-170 output and 
640x480 pixel resolution. Alternatively, camera 408 can comprise a CCD camera having 
CCIR output and 756x567 pixel resolution. 

In a particular embodiment of the invention, an infra-red illuminator 602 ("IR 
illuminator") is co-located with camera 108. IR illuminator 602 produces one or more 
beams of infrared light that is directed in the same direction as camera 108. IR illuminator 
602 comprises a surface that is ringed around the lens 606 of camera body 608, The surface 
of IR illuminator 602 contains a plurality of individual LED elements 604 for producing 
infrared light. The LED elements 604 are ananged in a spiral pattern on the IR illuminator 
602, Infrared filters that may be part of the camera 108 are removed or disabled to increase 
the camera's sensitivity to infrared light. 

According to an embodiment, digital video recordings of the patient in a session can 
be recorded via camera 108. The same camera 108 used for tracking patient movement can 
be used to record video unages of the patient for ftiture reference. A normal ambient light 
image sequence of the patient can be obtained in synchronization with the measured 
movement signals of markers 114. 

Figs. 7a and 7b depict an embodiment of a retro-reflective marker 700 that can be 
employed within the present invention. Retro-reflective marker 700 comprises a raised 
reflective surface 702 for reflecting light. Raised reflective surface 702 comprises a semi- 
spherical shape such that light can be reflected regardless of the input angle of the light 
source. A flat surface 704 surrounds the raised reflective surface 702. The underside of flat 
surface 704 provides a mounting area to attach retro-reflective marker 700to particular 
locations on a patient's body. According to an embodiment, retro-reflective marker 700 is 
comprised of a retro-reflective material 3M#7610WS available from 3M Corporation. In an 
embodiment, marker 700 has a diameter of approximately .5 cm and a height of the highest 
pomt of raised reflective surface 702 of approximately .1 cm. 

Fig. 8 depicts an q)paratus 802 that can be employed to manufacture retro-reflective 
markers 700. Apparatus 802 comprises a base portion 804 having an elastic ring 806 
affixed thereto. Elastic ring 806 is attached to bottom mold piece 808 having a bulge 
protruding from its center. A control lever 810 can be operated to move top portion 812 
along support rods 814. Top portion 812 comprises a spring-loaded top mold piece 814. 
Top mold piece 814 is formed with a semi-spherical cavity on its underside. In operation, a 
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piece of retro-reflective material is placed on bottom mold piece 808. Control lever 810 is 
operated to move top portion 812 towards base portion 804. The retro-reflective material is 
compressed and shaped between the bottom mold piece 808 and the top mold piece 814. 
The top mold piece 814 forms the upper exterior of the retro-reflective material into a semi- 
spherical shape. 

In an alternate embodiment, marker 1 14 comprises a marker block having one or 
more reference locations on its surface. Each reference location on the marker block 
preferably comprises a retro-reflective or reflective material that is detectable by an optical 
imaging apparatus, such as camera 108. 

Fig. 1 1 depicts an embodiment of a marker block 1 100 having a cylindrical shape 
with multiple reference locations comprised of retro-reflective elements 1 102 located on its 
surface. Marker block 1100 can be formed as a rigid block (e.g., from Styrofoam). Blocks 
made in this fashion can be reused a plurality of times, even with multiple patients. The 
retro-reflective elements 1 102 can be formed from the same material used to construct retro- 
reflective markers 1 14 of Figs. 7a and 7b. The marker block is preferably formed from a 
material that is light-weight enough not to interfere with normal breathing by the patient. 

A marker block can be formed into any shape or size, as long as the size, spacing, 
and positioning of the reference locations are configured such that a camera or other optical 
imaging apparatus can view and generate an image that accurately shows the positioning of 
the marker block. For example, Fig. 10 depicts an alternate marker block 1000 having a 
hemispherical shape comprised of a plurality of retro-reflective elements 1002 attached to 
its surface. 

The marker block can be formed with shapes to fit particular body parts. For 
example, molds or casts that match to specific locations on the body can be employed as 
marker blocks. Marker blocks shaped to fit certain areas of the body facilitate the 
repeatable placement of the marker blocks at particular locations on the patient. 
Alternatively, the marker blocks can be formed to fit certain fixtures that are attached to a 
patient's body. For example, a marker block can be formed within indentations and grooves 
that allow it to be attached to eyeglasses. In yet another embodiment, the fixtures are 
formed with integral marker block(s) having reflective or retro-reflective markers on them. 

An alternate embodiment of the marker block comprises only a single reference 
location/teflective element on its surface. This embodiment of the marker block is used in 
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place of the retro-reflective marker 406 to detect particular locations on a patient's body 
with an optical imaging apparatus. 

Further details regarding camera 108, markers 1 14, marker blocks, or procedures to 
optically measure physiological movements are described in copending U.S. Patent 
5 Application Ser, No. [Not Yet Assigned] (Varian Docket No. 98-34, Attorney Docket No. 
236/223) and U.S. Patent Application Ser. No. [Not Yet Assigned] (Varian Docket No. 98- 
33, Attorney Docket No. 236/224), filed concurrently herewith, both of which are hereby 
incorporated by reference in their entirety. 

Figure 12 is a block diagram that illustrates an embodiment of a computer system 
10 1900 upon which an embodiment of the invention may be implemented. Computer system 
1900 includes a bus 1902 or other communication mechanism for communicating 
information, and a processor 1904 coupled with bus 1902 for processing information. 
Computer system 1900 also includes a main memory 1906, such as a random access 
memory (RAM) or other dynamic storage device, coupled to bus 1902 for storing 
1 5 information and instructions to be executed by processor 1 904. Main memory 1 906 also 
may be used for storing temporary variables or other intermediate information during 
execution of instructions to be executed by processor 1904. Computer system 1900 further 
includes a read only memory (ROM) 1908 or other static storage device coupled to bus 
1902 for storing static information and instructions for processor 1904. A data storage 
20 device 1910, such as a magnetic disk or optical disk, is provided and coupled to bus 1902 
for storing information and instructions. 

Computer system 1900 may be coupled via bus 1902 to a display 1912, such as a 
cathode ray tube (CRT), for displaying information to a user. An input device 1914, 
including alphanumeric and other keys, is coupled to bus 1902 for communicating 
25 infonnation and command selections to processor 1904. Another type of user input device 
is cursor control 1916, such as a mouse, a trackball, or cursor direction keys for 
commxmicating direction information and command selections to processor 1904 and for 
controlling cursor movement on display 1912, This input device typically has two degrees 
of fi-eedom in two axes, a first axis (e.g., x) and a second axis (e.g., y), that allows the 
30 device to specify positions in a plane. 

The invention is related to the use of computer system 1900 for detecting and 
predictively estimating physiological cycles. According to one embodiment of the 
invention, such use is provided by computer system 1900 m response to processor 1904 
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executing one or more sequences of one or more instructions contained in main memory 
1906. Such instructions may be read into main memory 1906 from another computer- 
readable medium, such as storage device 1910. Execution of the sequences of instructions 
contained in main memory 1906 causes processor 1904 to perform the process steps 
described herein. One or more processors in a multi-processing arrangement may also be 
employed to execute the sequences of instructions contained in main memory 1906. In 
alternative embodiments, hard-wired circuitry may be used in place of or in combination 
with software instructions to implement the invention. Thus, embodiments of the invention 
are not limited to any specific combination of hardware circuitry and software. 

The term "computer-readable medium" as used herein refers to any medium that 
participates in providing instructions to processor 1904 for execution. Such a medium may 
take many forms, including but not limited to, non-volatile media, volatile media, and 
transmission media. Non-volatile media includes, for example, optical or magnetic disks, 
such as storage device 1910. Volatile media includes dynamic memory, such as main 
memory 1906- Transmission media includes coaxial cables, copper wire and fiber optics, 
including the wires that comprise bus 1902. Transmission media can also take the form of 
acoustic or light waves, such as those generated during radio wave and infrared data 
communications. 

Common forms of computer-readable media include, for example, a floppy disk, a 
flexible disk, hard disk, magnetic tape, or any other magnetic medium, a CD-ROM, any 
other optical medium, punch cards, paper tape, any other physical medium with patterns of 
holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any other memory chip or 
cartridge, a carrier wave as described hereinafter, or any other medium from which a 
computer can read. 

Various forms of computer readable media may be involved in carrying one or more 
sequences of one or more instructions to processor 1904 for execution. For example, the 
instructions may initially be carried on a magnetic disk of a remote computer. The remote 
computer can load the instructions into its dynamic memory and send the instructions over a 
telephone line using a modem. A modem local to computer system 1900 can receive the 
data on the telephone line and use an infrared transmitter to convert the data to an infrared 
signal. An infixed detector coupled to bus 1902 can receive the data carried m the infitired 
signal and place the data on bus 1902. Bus 1902 carries the data to main memory 1906, 
from which processor 1904 retrieves and executes the instructions. The instructions 
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received by main memory 1906 may optionally be stored on storage device 1910 either 
before or after execution by processor 1904. 

Computer system 1900 also includes a conmiunication interface 1918 coupled to bus 
1902. Communication interface 1918 provides a two-way data communication coupling to 
a network link 1920 that is connected to a local network 1922. For example, 
communication interface 1918 may be an integrated services digital network (ISDN) card or 
a modem to provide a data communication connection to a corresponding type of telephone 
line. As another example, communication interface 1918 may be a local area network 
(LAN) card to provide a data communication connection to a compatible LAN. Wireless 
links may also be implemented. In any such implementation, communication interface 
1918 sends and receives electrical, electromagnetic or optical signals that carry data streams 
representing various types of information. 

Network link 1 920 typically provides data communication through one or more 
networks to other devices. For example, network link 1920 may provide a connection 
through local network 1922 to a host computer 1924 or to medical equipment 1926 such as 
a radiation beam source or a switch operatively coupled to a radiation beam source. The 
data streams transported over network link 1920 can comprise electrical, electromagnetic or 
optical signals. The signals through the various networks and the signals on network link 
1920 and through communication interface 1918, which cany data to and from computer 
system 1900, are exemplary forms of carrier waves transporting the information. Computer 
system 1900 can send messages and receive data, including program code, through the 
network(s), network link 1920 and communication interface 1918. 

In the foregoing specification, the invention has been described with reference to 
specific embodiments thereof. It will, however, be evident that various modifications and 
changes may be made thereto without departing from the broader spirit and scope of the 
invention. For example, the operations performed by computer 1 10 can be performed by 
any combination of hardware and software within the scope of the invention, and should not 
be limited to particular embodiments comprising just a particular definition of "computer". 
The specification and drawings are, accordingly, to be regarded in an illustrative rather than 
restrictive sense. 
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